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SUMMARY 


In  calibrating  modern  inertial  quality  accelerometers  which  are  accurate 
to  a  few  parts  in  a  million,  extreme  care  must  be  exercised  in  measurement 
technique  and  interpretation  of  the  results.  This  Report  analyses  some  of  the 
tests  which  are  performed  currently  and  discusses  the  effects  of  errors  in  the 
test  equipment. 
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1  INTRODUCTION 

Although  an  assessment  laboratory  for  testing  inertial  navigation 

accelerometers  has  been  in  existence  for  over  20  years  at  RAE,  there  has  been  very 

little  documentation  of  testing  methods  and  philosophy.  The  major  work  entitled 

"Performance  characteristics  and  methods  of  testing  of  force-feedback 

accelerometers"  appeared  in  1967  and  while  this  gives  a  comprehensive  account  of 

the  underlying  physical  principles,  the  test  methods  described  largely  relate  to 

the  equipment  available  at  the  time,  with  emphasis  on  manual  recording  of  data  and 

.  4  .  . 

an  accuracy  of  the  order  of  1  part  in  10  .  With  the  advent  of  digital  data 
acquisition  systems  and  rapid  machine  processing  of  data,  many  more  tests  are  now 
performed  than  hitherto  possible  and  some  existing  test  methods  have  been 
improved.  Furthermore,  accelerometers  are  now  available  from  a  number  of 
different  manufacturers  who  claim  accuracies  in  the  order  of  a  few  parts  in  10^. 

This  Report  is  intended  to  provide  a  description  and  analysis  of  some  of  the 
test  methods  currently  employed  at  RAE  and  to  indicate  where  possible  the 
magnitudes  and  sources  of  error  likely  to  be  encountered.  Little  will  be  said 
about  the  operating  principles  of  force  feedback  accelerometers  or  the  character¬ 
istics  of  the  servo  electronics  since  these  aspects  have  been  discussed  in  detail 
1  2 

elsewhere  ’  .  However,  it  is  assumed  that  the  instruments  are  pendulous  devices 
with  a  single  degree  of  freedom,  in  which  the  proof  mass  is  captured  in  a 
sufficiently  tight  servo  loop  so  that  errors  associated  with  transverse 
accelerations  are  small.  The  pendulum  is  supported  by  means  of  flexure  hinges  or 
pivot  and  jewel  arrangements.  Non-pendulous  devices  or  those  with  more  than  one 
degree  of  freedom  are  not  so  commonly  available,  although  much  of  the  present 
analysis  may  be  extended  to  include  these. 

There  are  several  tests,  for  example  those  which  measure  the  properties  of 
the  flexure  or  pick-off,  which  indicate  that  the  instrument  has  been  built  to  the 
required  design  standard  and  is  capable  of  functioning  correctly.  They  do  not, 
however,  measure  performance  under  normal  operating  conditions  and  are  open  loop 
tests  requiring  access  to  the  pick-off.  Many  of  today's  accelerometers  are 
supplied  as  packaged  units  comprising  integrated  electronics  and  mechanical  sensor 
and  no  provision  is  made  for  the  user  to  obtain  access  to  monitor  test  points 
within  the  unit.  Since  this  Report  is  concerned  mainly  with  operational  perform¬ 
ance  and  such  tests  have  already  been  described  adequately  in  Ref  1 ,  they  will  not 
be  considered  further. 

In  most  analogue  instruments,  the  output  quantity  is  the  current  flowing  in 
the  restoring  coil  which  is  usually  measured  directly.  Some  accelerometers. 
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however,  are  designed  to  give  a  voltage  output  and  this  may  be  influenced  by  the 
stability  of  the  built-in  load  resistor.  Digital  instruments  also  exist  which 
employ  a  pulse  torque  rebalance  loop  and  the  output  is  then  a  series  of  pulses. 

So  far,  no  instruments  of  this  type  have  been  tested  in  this  laboratory  although 
most  of  the  analysis  which  follows  is  still  applicable.  The  output  is  assumed  to 
be  some  reasonable  function  of  the  input  acceleration  and  will  almost  certainly 
be  affected  by  changes  in  the  environmental  or  operating  conditions. 

2  PRINCIPLES  OF  TESTING 

It  is  standard  practice  to  perform  a  large  part  of  accelerometer  assessment 
using  earth's  gravity  as  the  input.  In  spite  of  the  fact  that  the  range  of  most 
instruments  (typically  ±20g)  is  much  greater  than  the  ±lg  available  from  earth's 
field  testing,  a  large  amount  of  useful  information  can  be  gained.  Experience  in 
several  test  laboratories,  mainly  in  the  USA,  has  indicated  that  extrapolated 
results  from  ±lg  tests  usually  compare  favourably  with  those  using  a  precision 
centrifuge  to  provide  the  necessary  high-g  levels.  In  any  case,  it  is  very 
difficult  to  eliminate  or  correct  for  all  the  errors  associated  with  practical 
centrifuges. 

In  order  to  evaluate  and  compare  units  from  the  present  generation  of 

precision  accelerometers,  both  output  current  and  input  acceleration  must  be 

known  to  a  high  degree  of  accuracy.  The  instrument  generally  used  at  RAE  to 
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measure  output  is  a  digital  integrating  milliammeter  whose  resolution  can  be  as 
good  as  the  equivalent  of  0.1  yg  for  sufficiently  long  integration  periods.  Such 
accuracy  is  rarely  required  except  for  the  most  precise  stability  tests,  and  in 
most  cases  resolutions  of  1  or  10  yg  are  accepted  so  that  shorter  integration 
times  can  be  employed.  To  determine  the  input,  it  is  usual  to  place  the  accelero¬ 
meter  at  precisely  known  angles  with  respect  to  the  gravity  vector  by  mounting  it 
securely  on  a  fixture  attached  to  the  faceplate  of  an  accurate  dividing  head. 

The  Leitz  dividing  heads  used  at  RAE  and  several  other  test  centres  have  an 
angular  setting  accuracy  of  about  ±1  arc  second  which  corresponds  to  about  5  yg 
when  considered  as  a  deviation  from  the  horizontal  plane.  The  precision  with 
which  the  dividing  head  is  initially  set  up  and  the  orthogonality  of  the  mounting 
fixture  are  important  considerations  and  the  errors  associated  with  these  are 
discussed  in  section  4. 

Account  must  be  taken  of  any  external  disturbances  which  are  likely  to 
affect  the  performance  of  a  test.  Temperature  is  usually  the  dominant  influence 
and  since  it  is  difficult  to  control  the  test  environment  to  better  than  1°C,  the 
temperature  at  the  accelerometer  should  be  monitored  continuously  and  used  to 


correct  the  raw  data  if  necessary.  In  this  laboratory,  nickel  resistance 
temperature  sensors  are  generally  employed  in  conjunction  with  digital  thermometers 
having  a  resolution  of  0.01°C  in  the  range  -10  to  +  120°C^  .  It  is  usually  found 
that  there  is  excellent  correlation  between  temperature  and  accelerometer  output  at 
a  fixed  attitude  and  thus  corrections  can  be  made  with  a  high  degree  of  confidence. 

Tests  should  be  repeated  a  number  of  times  under  the  same  conditions  to 
enable  an  estimate  to  be  made  of  the  test  uncertainty.  An  e*  'mple  is  the  four- 
point  test  in  which  the  output  is  measured  at  the  two  Og  positions  and  the  two  Ig 
positions  as  the  accelerometer  is  rotated  about  one  of  its  non-sensitive  axes  (see 
section  3.2).  It  is  usual  practice  to  perform  this  test  four  vimes  with  a  repeat 
reading  of  the  +lg  output  during  each  revolution.  The  mean  values  and  standard 
deviations  of  the  parameters  to  be  measured  are  calculated  together  with  the  in¬ 
test  rms  deviation  which  estimates  the  test  accuracy.  Should  this  value  be 
unacceptably  large,  then  the  results  should  be  rejected  if  there  is  any  suspicion 
that  the  operator  or  measuring  equipment  were  at  fault.  It  is  damaging  to  the 
reputation  of  the  laboratory  and  unfair  to  the  manufacturer  if  his  accelerometer 
is  blamed  for  imprecise  or  faulty  test  methods. 

It  is  also  good  practice  to  determine  performance  parameters  using  different 
test  procedures.  For  example,  temperature  coefficients  can  be  measured  by  perform¬ 
ing  tests  at  fixed  temperatures  within  a  range,  or  by  establishing  steady  conditions 
at  a  high  temperature  and  recording  the  cooling  curves  at  different  attitudes. 
Similar  results  give  added  confidence  in  the  values  obtained. 

The  introduction  of  automatic  data  acquisition  and  processing  has  resulted  in 
a  number  of  significant  improvements  in  accelerometer  assessment,  which  include 

(a)  less  reliance  on  operator  skill  and  virtual  elimination  of  arithmetical 
mistakes ; 

(b)  a  greatly  increased  number  of  tests  and  hence  better  estimates  of  scatter  and 
drift; 

(c)  more  reliable  output  stability  information  over  long  periods  of  time  with 
automatic  correction  for  temperature; 

(d)  the  ability  to  analyse  test  data  obtained  at  a  large  number  of  intermediate 
g-levels  and  hence  more  accurate  measurements; 

(e)  the  collection  of  statistical  information  which  enables  levels  of  confidence 
to  be  established  and  predictions  to  be  made  about  long  term  behaviour; 


(f)  automatic  graph  plotting. 
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At  present,  the  output  from  digital  instruments  (milliammeters ,  thermometers, 
etc)  is  in  parallel  BCD  format  and  the  data  from  a  number  of  instruments  is 
collected,  reformatted  and  then  transferred  on  to  punched  paper  tape  for  subse¬ 
quent  computer  analysis.  Further  improvements  are  envisaged  in  which  the  instru¬ 
mentation  would  be  linked  together  via  a  standard  interface  bus  to  a  digital 
controller  and  the  paper  tape  would  be  replaced  by  some  form  of  magnetic  storage, 
either  cartridge  or  disc.  A  more  distant  aim  is  to  automate  the  dividing  head 
setting  to  enabli  the  majority  of  testing  to  be  performed  without  manual 
intervention. 

3  DERIVATION  OF  ACCELEROMETER  PARAMETERS 
3. 1  The  model  equation 

The  output  of  an  ideal  accelerometer,  expressed  in  units  of  '  g'  can  be 
stated  as 


I  =  Vs  (3-1) 

where  Kj  is  the  scale  factor  and  ag  is  the  component  of  acceleration  along  the 
sensitive  axis.  In  practice,  the  situation  is  more  complex.  There  will  be,  in 
general,  a  non-zero  output  when  no  accelerations  are  acting  (bias);  small  misalign¬ 
ments  between  the  nominal  and  true  accelerometer  axes  will  be  present;  the  output 
may  not  be  a  strictly  linear  function  of  acceleration  and  there  may  be  s.  slight 
sensitivity  to  accelerations  perpendicular  to  the  sensitive  axis.  There  may  also 
be  additional  errors  due  to  cross  coupling.  The  choice  of  model  equation  depends 
to  a  certain  extent  on  the  type  of  accelerometer  to  be  tested.  Digital  instruments 
employing  pulse  rebalance  electronics  sometimes  have  asymmetrical  scale  factors  and 
biases  because  of  the  separate  circuits  used  to  generate  the  positive  and  negative 
pulses.  The  same  may  be  true  of  accelerometers  which  are  not  symmetrical  about 
their  pendulum  axis,  for  example  those  with  only  one  restoring  coil  and  magnet 
rather  than  two.  However  an  equation  which  has  been  found  to  model  realistically 
the  output  of  most  accelerometers  tested  at  RAE  is 


I  =  Kn  +  K  a  +  K0a^  +  K,a^  +  6  a  +  +  K  a  a  +  K 

0  Is  2  s  3  s  1  p  2  n  sp  s  p 


shas3h 


(3-2) 


where  is  the  bias,  K2  and  are  non-linearity  coefficients,  6^  and 

<$2  are  transverse  sensitivity  coefficients,  and  Kg^,  Kgh  are  the  cross-coupling 
coefficients.  The  components  a  ,  a  and  a,  are  the  accelerations  along  the 

s  p  tl 

true  sensitive,  pendulum  and  hinge  axes  respectively  according  to  an  orthogonal 
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system  in  which  s  *  g  =  h  (see  Fig  la).  The  coat j.ic' ents  in  equation  (3~2)  may 
themselves  be  functions  of  the  environment  or  operating  conditions.  They  are 
therefore  likely  to  change  with  changes  in  temperature,  amplifier  gain,  mechanical 
shock,  etc. 

In  order  to  determine  the  coefficients  in  equation  (3-2) ,  the  accelerometer 
is  set  in  various  orientations  with  respect  to  earth's  gravity  and  the  output  is 
measured  under  constant  environmental  conditions.  This  orientation  is  determined 
with  respect  to  an  orthogonal  set  of  reference  (or  nominal)  axes  which  are  usually 
defined  by  the  accelerometer  structure.  Most  accelerometers  provided  with  a 
lapped  mounting  surface  containing  three  or  four  fixing  holes,  two  of  which  are 
nominally  parallel  to  the  hinge  or  pendulum  axes.  In  a  typical  case,  we  may 
define  the  nominal  sensitive  axis  as  being  perpendicular  to  the  plane  of  the 
mounting  surface  and  the  nominal  hinge  axis  as  being  in  this  plane  and  joining  the 
centres  of  the  two  appropriate  fixing  holes.  The  nominal  pendulum  axis  is  then 
perpendicular  to  both  of  these.  We  shall  use  primed  quantities  (a^,  a^,  a^)  to 
denote  acceleration  components  referred  to  the  nominal  coordinate  frame. 

Consider  an  orthogonal  set  of  axes  fixed  in  space  (Fig  lb)  such  that  y^ 
is  vertically  down  (parallel  to  the  gravity  vector),  z^  is  horizontal  and  in  the 
vertical  plane  containing  the  spindle  of  the  dividing  head,  and  x^  is  perpendic¬ 
ular  to  y„  and  z.  .  Then  a  =  a  =0  and  a  =  1  .  For  the  time  being  ye 
0  0  *0  z0  *0 
shall  make  four  assumptions,  the  validity  of  which  will  be  tested  in  the  next 


section: 


the  dividing  head  is  error-free; 

its  axis  of  rotation  is  accurately  level; 


(c)  the  accelerometer  is  mounted  so  that  one  of  its  nominal  axes  is  parallel  to 
the  axis  of  rotation  of  the  dividing  head; 

(d)  in  the  reference  (0  =  0)  position,  the  other  nominal  axes  are  precisely 
horizontal  and  vertical. 

Then,  if  x',  y'  and  z'  are  the  nominal  accelerometer  axes,  and  the  dividing 
head  is  rotated  through  an  angle  0  ,  then 


cos  0 


sin  0 


-  sin  0  00 


cos  0  01 


(3t3) 


10 


t&MiVcir 
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Note  that  the  dividing  head  is  mechanised  so  that  the  direction  of  increasing  0 
is  as  shown  in  Fig  1 ,  and  that  accelerometers  are  conventionally  tested  by 
rotating  them  in  the  sense  +lg  through  the  hinge-up  and  -lg  positions  to  the 
hinge-down  position.  If  the  accelerometer  is  mounted  so  that  its  nominal  hinge 
axis  is  parallel  to  the  axis  of  rotation,  which  is  the  most  common  configuration 
for  testing  and  referred  to  here  as  mounting  position  i ,  then 


a'  = 

a*  =  -  sin  0 

(3-4) 

X 

s 

a*  = 

y 

a'  =  cos  0 

P 

(3-5) 

a'  * 
z 

a-  =  0  . 

(3-6) 

In  general,  the  nominal  and  true  axes  will  not  be  coincident  and  an  ortho¬ 
gonal  transformation  can  be  set  up  which  relates  the  two  frames  of  reference. 

Let  the  transformation  from  true  to  nominal  axes  be  described  by  three  sequential 
rotations;  the  first  through  an  angle  a  about  the  (true)  hinge  axis;  the  second 
through  an  angle  0  about  the  pendulum  axis  generated  by  the  first  rotation  and 
the  last  through  an  angle  y  about  the  now  nominal  sensitive  axis.  This  yields 
the  transformation  matrix,  [C]  = 

cos  a  cos  0 

-  sin  a  cos  y  +  cos  a  sin  0  sin  y 
sin  a  sin  y  +  cos  a  sin  0  cos  y 
in  the  equation 

a'  =  [c]  a 

where  a  and  a’  are  vectors  representing  the  true  and  nominal  axes  respectively. 
Since  the  matrix  is  orthogonal,  the  inverse  transformation  is 

a  =  [Cf1  a'  =  [C  ]T  a'  (3-7) 

where  [c]  is  the  transposed  matrix.  The  transformation  has  been  carried  out 
in  this  way  so  that  the  generally  accepted  misalignment  angles  are  introduced  in 
their  correct  sequence;  a  and  0  are  the  misalignments  of  the  true  and  nominal 
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sensitive  axes  about  the  hinge  and  pendulum  axes  respectively  and  y  represents 
a  misalignment  about  the  sensitive  axis.  For  an  accelerometer  to  be  of  inertial 
quality  the  misalignments  must  be  small.  If  the  angles  are  less  than  about  3  arc 
minutes,  then  an  accuracy  of  1  part  in  109  is  retained  if  we  take  first  order 
approximations .  Hence 


1  ~  a  3 

a  1  -  y 


Y 


1 


(3-8) 


and  using  equations  (3  —4)  and  (3~8),  we  have 

a  =  -  sin  0  -  a  cos  0  (3~9) 

s 

a  =  -  a  sin  0  f  cos  0  (3“10) 

P 

a^  =  8  sin  0  +  y  cos  0  .  (3~11) 


Substituting  equations  (3~9)  to  (3~11)  into  the  model  equation  (3-2)  and 
neglecting  products  of  small  quantities  such  as  K^a,  K^a,  6a,  etc;  we  have 
(where  the  additional  subscript  is  used  to  denote  mounting  position  1) 

I01  =  K01  ”  K11  S^n  9  +  K21  s^n2e  “  K3i  sin39 

+  (6  -  K..a)  cos  0  -  (K  sin  20  .  (3-12) 

ill  sp 

Thus  small  misalignments  about  axes  other  than  that  perpendicular  to  the  plane  of 
rotation  do  not  produce  any  significant  errors  in  the  measurement  of  the 
coefficients . 

A  similar  relationship  holds  for  mounting  position  2  in  which  the  nominal 
pendulum  axis  is  maintained  parallel  to  the  axis  of  rotation.  In  this  case 


a'  = 

a:  = 

-  sin  0 

X 

h 

a'  - 

a’  = 

cos  0 

y 

s 

a’  = 

a'  = 

0 

z 

P 
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and  the  components  of  acceleration  along  the  true  axes  are. 


a  =  cos  0-8  sin  0 
s 


a  =  a  cos  9  +  y  sin  0 
P 


=  -  8  cos  0  -  sin  0 


Equation  (3.2)  becomes 


I02  =  *02  +  *12  C0S  6  +  *22  cos2e  +  K32  cos3e 


-  (6^  +  sin  *■’  ~  sin  26 


(3-13) 


where  the  extra  subscript  is  used  to  den  ite  mounting  position  2.  A  third  ortho¬ 
gonal  mounting  position  is  possible,  ie  that  in  which  the  sensitive  axis  is 
parallel  to  the  axis  of  rotation.  However,  this  arrangement  is  not  often  used  in 
practice  since  the  results  are  much  more  sensitive  to  errors  in  the  dividing  head 
than  they  are  in  the  two  standard  mounting  positions  (see  section  4) .  Assuming  no 
errors,  the  output  in  position  3  is  given  by 


I03  =  KQ3  -  (6]  -  Kj3a)  sin  0  +  (62  +  K.^8)  cos  0  .  (3-14) 


3.2  Four-point  tests 

The  main  performance  parameters  can  be  deduced  from  simple  tests  in  which 
the  accelerometer  output  is  measured  at  each  of  four  cardinal  positions.  Assuming 
the  unit  is  in  mounting  position  1  and  using  HU  and  HD  to  denote  the  zero-g 
hinge-up  and  hinge-down  attitudes  respectively,  we  have  the  following  set  of 


equations. 

0g  (HU) 

xo 

=  K01 

+  (6j  -  Kj , a) 

-lg 

I90 

=  Ko, 

-K.l  +K21  " 

*31 

0g  (HD) 

o 

00 

M 

=  *01 

-  (5j  -  Kjjo) 

+  lg 

:270 

=  Koi 

+K11  *21  + 

*31 

11 


Hence 


K01  ^(I0  +  I180) 

=  Bl, 

(bias  at  Og) 

(3-15) 

K01 

+  K21  =  *(I90  +  I270) 

"  B21 

(bias  at  lg) 

(3-16) 

K21  =  B21  "  B11 

-  BD 

(bias  discrepancy) 

(3-17) 

6, 

JK 

II 

K#— 

/-\ 

M 

O 

1 

M 

CO 

O 

v-/ 

=  AE 

(axis  error) 

(3-18) 

K11 

+  K31  *(I90  "  I270) 

“  Si, 

(scale  factor)  . 

(3-19) 

There  are  two  points  to  note  here.  The  first  is  that  it  is  not  possible  to 
separate  the  effects  of  any  inherent  transverse  sensitivity  from  true  axis  mis¬ 
alignments.  What  we  measure,  therefore,  is  a  lumped  parameter  known  as  axis 
error,  or  more  correctly,  the  component  of  axis  error  in  the  plane  containing  the 
sensitive  and  pendulum  axes.  The  second  point  is  that  the  difference  between  the 
+lg  and  -lg  readings  gives  the  scale  factor  plus  a  small  additional  term  due  to 
third  order  non-linearity.  In  most  cases,  this  quantity  is  so  small  as  to  be 
negligible,  and  the  value  measured  represents  scale  factor  to  an  adequate  degree 
of  accuracy.  In  mounting  position  2,  the  relationships  are 


o 

to 

II 

*(I90  +  I270) 

-  31 2 

K02+K22  " 

*(IQ  +  I180) 

"  B22 

K22  " 

B22  "  B,2 

=  bd2 

62  ~  K12B  = 

^(I90  ”  1270) 

=  HAE 

K12+K32  " 

N*— 

/-s 

M 

O 

1 

M 

5 

o 

sf2 

(bias  at  Og) 

(3-20) 

(bias  at  lg) 

(3-21) 

(bias  discrepancy) 

(3-22) 

(hinge  axis  error) 

(3-23) 

(scale  factor) 

(3-24) 

Again,  equation  (3~23)  shows  that  what  we  measure  is  a  parameter  described  as  the 
component  of  axis  error  in  the  plane  containing  the  sensitive  and  hinge  axes. 

Some  authorities  refer  to  this  as  the  hinge  axis  error. 

The  values  of  the  other  parameters  should  be  the  same  whether  measured  in 
mounting  position0  1  or  2.  In  practice,  because  of  dividing  head  errors  and 
149  other  small  misalx^  .ments  considered  negligible  in  the  above  analysis,  small 
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differences  may  be  found,  but  should  be  insignificant  in  any  reasonable 
application.  It  should  also  be  noted  that  the  cross-coupling  coefficients  do  not 
appear  in  the  analysis  of  four-point  tests. 

3 . 3  Mnltiposition  tests 

The  four-point  test  only  provides  information  at  the  zero-g  and  two  extreme 
lg  levels.  While  this  is  not  usually  a  serious  limitation  for  routine  day-to-day 
assessment,  better  estimates  of  the  performance  parameters  are  obtained  if  they 
are  derived  from  data  taken  at  a  large  number  of  intermediate  g-levels  as  well  as 
the  extreme  values.  In  a  multiposition  test,  the  output  is  recorded  at  different 
angles  and  regression  analysis  is  used  to  fic  the  data  to  the  chosen  model 
equation.  An  additional  advantage  is  that  the  validity  of  the  model  can  be  tested 
by  examining  a  plot  of  the  residuals,  and  if  necessary  the  data  can  be  fitted  to  a 
different  model,  for  example  one  using  asymmetrical  scale  factors. 

In  the  following  analysis  we  shall  assume  that  the  model  equation  (3-2)  is 
appropriate  and  that  the  accelerometer  is  in  mounting  position  1 .  Then  if  1^ 
denotes  the  best  fit  value  of  output  current  determined  from  the  theoretical  curve 
at  position  k  ,  equation  (3-12)  becomes 

2  3 

Ik  *  Kqj  -  Kj j (~  sin  0^+  a  cos  0^)  +  K2]  sin  0^  -  K^j  sin  0fc 

+  6,  cos  0,  -  JK  sin  20,  .  (3-25) 

1  k  sp  k 

If  I  denotes  the  measured  output  at  position  k  ,  then  we  use  the  method  of 

tv 

least  mean  squares  to  minimise  the  sum  of  the  squares  of  the  residuals  with 
respect  to  each  unknown  coefficient  we  wish  to  determine;  %e  we  need  to  minimise 

I  <  -  £«*  -  v2 

all  k 


m 


+  K, .  (sin  0,  +  a  cos  0,  )  - 
11  k  k 


(3-26) 


with  respect  to  ,  Kj  j ,  K^j ,  K^j ,  <S,  K  and  a  .  We  therefore  form  the  partial 
derivatives  of  equation  (3-26)  with  respect  to  each  unknown  and  set  these  equal  to 
zero.  It  might  be  argued,  with  some  justification,  that  the  0^  should  be 
chosen  to  produce  approximately  equal  increments  of  applied  acceleration.  If 
this  scheme  is  employed,  however,  the  computations  are  extremely  involved, 
requiring  the  inversion  of  an  m  *  m  matrix  where  m  is  the  number  of  unknown 
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coefficients.  To  simplify  the  data  reduction  procedure,  we  shall  assume  that 

0^  =  kir/2n  where  n=  1,  2,  3,  .  and  k  =  0,  1,  2  .  (4n  -  1).  This 

means  that  we  perform  an  N-position  test,  where  N  is  a  multiple  of  4,  and  we 
use  equal  angular  increments. 

Differentiating  equation  (3-26)  with  respect  to  Kqj,  we  have 


where  all  summations  are  taken  from  k  =  0  to  k  =  4n  -  1  . 

Therefore  "  Kl  1  Z  ^  ‘  “K"  Z  ^ 

+  K21  Z  Si"2©  "  K3>  Z  Si"3® 

+  si  Z-  ®  ■  !ksp  Z sin  ^  • 

This  simplifies  to 

Zk  “  4nKoi  *  2l,K2i  •  <3’27) 

When  we  differentiate  equation  (3-26)  with  respect  to  all  the  other  unknown 
quantities  and  simplify,  we  obtain  the  following  expressions; 

Z1*  sin  ©  +  “  Z  **  c°s  ©  '  ■  2nK»  "  T  K3I  <3-28) 

Z\-n2©  "  2nK0.  +TK21  (3-29) 

Zk  Sin3(^  =  'f  Kl!  --TK3,  «-3« 
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«  ©  ■  2"<5,  -Kn»> 

V'  I,  sin  (— )  =  -  nK 

/  ,  k  4'  sp 


(3-31) 


(3-32) 


These  expressions,  with  the  exception  of  equation  (3-31),  are  not  valid  for 
n  =  1  ,  ie  a  four-point  test.  Note  that  equation  (3-31)  is  obtained  when  we 
differentiate  with  respect  to  either  6^  or  a  ,  again  showing  that  it  is  not 
possible  to  separate  these  coefficients.  Equations  (3-28)  to  (3-32)  reduce  to 


(3-33) 

(3-34) 

(3-35) 

13-36) 

(3-37) 

(3-38) 


The  above  equations  represent  the  best  estimates  of  the  coefficients  in  equation 

(3-12)  when  measured  from  an  N-position  test,  where  N  =  8,  12,  16,  20,  ....,  . 

Since  K_.  and  K  should  be  very  small,  it  is  unlikely  that  any  test  over 
j  *  sp 

the  narrow  ±lg  range  would  produce  meaningful  values.  Should  they  in  fact  be 
significantly  greater  than  the  test  uncertainty  as  measured  by  the  standard 
deviation  of  the  residuals,  then  the  test  would  almost  certainly  indicate  a  poor 
or  faulty  accelerometer. 

An  identical  analysis  holds  for  mounting  position  2  which  allows  measure¬ 
ments  to  be  made  of  the  hinge  axis  error  and  of  the  sensitive/hinge  axes  cross 
coupling  coefficient. 
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4  DIVIDING  HEAD  AND  LEVELLING  ERRORS 


It  has  been  assumed  so  far  that  there  have  been  no  errors  associated  with 
the  dividing  head  and  mounting  fixture.  In  practice,  of  cou-se,  this  is  not  the 
case  and  in  this  section,  an  attempt  is  made  to  identify  those  errors  which  exist 
and  to  determine  the  magnitudes  that  can  be  tolerated  while  still  maintaining  the 
required  measurement  accuracy. 

Let  us  assume  that,  instead  of  lying  along  the  z^  (horizontal)  axis,  the 
axis  of  rotation  of  the  dividing  head  is  inclined  at  a  small  angle  to  the 

horizontal  (see  Fig  2).  Then  the  acceleration  components  along  the  x^  y 
axes  are 


-  sm  e. 


o  rn 


sm  e.  a 

1  yr 


cos  s.  I  a 
1 1 


(4-1) 


where  3y^>  az  j  =  {0»  Of  .  If  the  0  =  0  position  is  not  precisely 

horizontal,  but  misaligned  by  a  small  angle  »  the  new  acceleration  components 
are  given  by 


0  a 


0  a 


1  I  a 

J  L  z' 


(4-2) 


This  is  the  situation  when  the  accelerometer  has  been  nominally  levelled  (6  =  0°) , 
where  ja^,  ay^’  az2}  are  the  acceleration  components  along  the  3  nominal  axes. 

As  the  dividing  head  spindle  is  rotated  through  any  angle  0  ,  the  components  are: 


cos  e 


sin  9 


-  sin  9 


cos  e 


(4-3' 
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and  combining  equations  (4-1),  (4-2)  and  (4-3)  we  have 


cos  (8  +  e^) 


sin  (6  +  e2) 


sin  (8  +  z^)  cos  g^ 

-  sin  (8  +  Zy)  sin 

a 

.  xo 

cos  (0  +  cos  ej 

cos  (8  +  tj)  sin  Gj 

a 

yo 

-  sin 

cos  G  j 

a 

_  zo_ 

Hence 

ar  =  -  sin  (0  +  g  )  cos  g 

X3  Z  I 

a  =  cos  (8  +  znj  cos  e, 
y3  2  1 

a  =  sin  e,  . 

Z3  1 


(A-4) 

(4-5) 

(4-6) 


This  is  usually  as  far  as  the  analysis  is  taken  and  it  shows  that,  provided 
the  accelerometer  is  mounted  in  a  position  such  that  its  nominal  sensitive  axis 
rotates  in  the  xy  plane,  only  the  misalignment  contributes  significantly  to 

errors  in  measuring  the  performance  parameters,  specifically  the  axis  error.  In 
practice,  care  is  taken  to  reduce  to  within  about  2  arc  seconds,  whereas  an 

angular  misalignment  of  a  few  arc  minutes  in  can  be  tolerated  before  notice¬ 

able  measurement  errors  are  introduced.  When  tha  sensitive  axis,  however,  is 
aligned  along  the  z-direction,  dividing  head  anc  mounting  block  inaccuracies  are 
significant  as  we  shall  see  below. 

Two  further  misalignments  can  be  introduced  as  follows.  An  accelerometer 
mounting  fixture  is  usually  designed  to  allow  one  of  the  nominal  axes  to  be 
aligned  parallel  to  the  axis  of  rotation  of  the  dividing  head.  As  a  result  of 
machining  limitations  or  setting  up  inaccuracies,  there  can  be  a  misalignment 
between  this  accelerometer  nominal  axis  and  the  dividing  head  spindle  about  the 
x-direction  and  an  additional  misalignment  g^  about  the  y-direction  (see  Fig  2) . 
The  way  in  which  these  errors  arise  depends  upon  the  detailed  fixing  arrangements, 
but  in  a  typical  case  a  plane  lapped  surface  is  provided  with  a  pattern  of  fixing 
holes  on  to  which  the  accelerometer  is  bolted.  A  skew  error  g^  will  be;  present 
if  this  plane  is  not  precisely  perpendicular  to  the  dividing  head  faceplate;  or 
indeed  if  the  faceplate  is  not  perpendicular  to  the  spindle  axis.  If,  in  addition, 
the  mounting  hole  configuration  does  not  permit  precise  alignment  between  the 
required  accelerometer  nominal  axis  and  the  axis  of  rotation,  then  a  skew  error 
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will  arise.  In  some  circumstances,  can  be  relatively  large  (~  20  arc 

minutes)  since  there  is  usually  a  certain  slackness  in  the  fixing  holes  and  the 
adjustment  is  done  by  eye.  It  is  worth  noting  also  that  significant  errors  may 
still  be  present  in  a  system  which  appears  to  be  perfectly  aligned.  For  example, 
the  accelerometer  mounting  surface  may  be  adjusted  to  be  precisely  horizontal  at 
one  particular  dividing  head  setting,  say  0=0°,  and  it  would  be  tempting  to 
believe  that  all  errors  had  been  reduced  to  zero.  However,  the  same  result  would 
have  been  obtained  if  the  axis  of  rotation  had  not  been  level  and  the  mounting 
surface  had  been  skewed  by  an  equal  but  opposite  amount,  ie  .  In  this 

case  gross  errors  might  result  in  the  measurements  of  the  performance  parameters. 
It  is  therefore  important  to  ensure  that  the  dividing  head  is  set  up  to  a 
sufficient  degree  of  accuracy  before  attaching  any  mounting  fixtures,  and  that 
reasonably  tight  tolerances  are  applied  in  the  manufacture  of  such  fixtures. 

If  x4>  y^»  z4  anc*  x^ ,  y,.,  z j.  are  the  axes  after  rotation  through 
and  ,  we  have  for  the  acceleration  components: 


and 


0  cos  sin 


0  -  sin  cos  £3 


cos  e. 


sin  e.  0 

4 


0  -  sin  e. 


cos  e. 


(4-7) 


(4-8) 


This  final  transformation  gives  the  components  of  acceleration  along  the 
accelerometer  nominal  axes  after  taking  into  account  all  the  levelling  and  mount¬ 
ing  misalignments  which  should  remain  constant  during  a  test  in  which  the 
instrument  is  not  removed  from  its  fixture.  No  account  has  been  taken  of  test 
plinth  instability  caused  by  random  vibrations  or  temperature  changes,  but  these 
effects  are  largely  unpredictable  and  should  be  negligibly  small  in  any  worthwhile 
test  facility.  It  has  also  been  assumed  that  the  dividing  head  is  of  such 
precision  that  backlash  in  the  gearing  can  be  ignored. 
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Since  primed  quantities  are  being  used  to  denote  accelerations  along  the 

nominal  axes,  we  have  <a  ,  a  ,  a  >=  <a',  a'»  a'  ?  and  equations  (4-1)  to 

y=  zc/  1  x  ”  z  J 


5  3  5  5 

(4-8)  can  be  combined  to  give 


a*  =  -  sin  (0  +  O  cos  e,  cos  e,  +  cos  (0  +  e„)  cos  e,  sin  e  sin  e 

x  2  14  2  1  o  4 


+  sm  e,  cos  £„  sm  e, 
1  3  4 


(4-9) 


a^  =  cos  (0  +  e^)  cos  cos  £^  ~  s^n  e]  s^n  £3  (4-10) 

a'  =  -  sin  (0  +  e.)  cos  e.  sin  e. 

Z  Z  1  4 


Note  that  if  there  are  no  errors,  these  reduce  to  equation  (3-3). 

To  retain  an  accuracy  of  1  part  in  10^,  we  may  take  first  order  angle 
approximations  only  if  the  angles  are  less  than  about  3  arc  minutes.  Although 
this  is  true  of  z ^  (which  we  always  take  care  t>  reduce  to  near-zero)  it  is  not 
necessarily  the  case  with  e^,  and  .  However,  these  angles  will  almost 
certainly  be  less  than  i  and  there  is  no  loss  of  accuracy  if  we  make  the 
approximations 


-  cos  (0  +  e^)  cos  sin  cos 

-  sin  e,  cos  £ n  cos  e,  .  (4-11) 

1  3  4 


sin  £ .  ~  £ . 
l  l 


cos  £ 


i  ■  1  -  <1 


2  (i  =  1,  3,  4) 


Equations  (4-9),  (4-10)  and  (4-11)  then  become 


2  2 


sin  0  -  (e^  ~  e3e^)  cos  ®  +  eie4 


/  2  2  v 


,»  =  - 


sin  0  + 


cos  0  - 


(4-12) 


(4-13) 


=  - 


£,  sin  0  —  £„  cos  0  -  £, 


(4-14) 
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To  determine  the  acceleration  components  along  the  true  axes,  we  use  the  trans¬ 
formation  (3-8).  Assuming  first  that  the  accelerometer  is  mounted  in  position  1, 
ie  its  hinge  axis  is  nominally  parallel  to  the  axis  of  rotation,  we  have,  after 
making  the  appropriate  small  angle  approximations: 


2  2 


=  “  ^  -  Y  ~  ~2  ^  sin  6  ~  (“ 


a  = 
s 


a  =  -  (a  +  e„l  sin  8  + 

p  v  2J 


+  e2  -  e.^)  cos  6 


+  e  e  (4-15) 

1  4 


/  e2  e2  \ 

'T't)cos  0  *  sIe3 


(4—1 6) 


3  =  fd  -  e.l  sin  0  +  fy  -  e.l  cos  0  -  e. 

n  v  4'  v  3'  1 


(4-17) 


Hence,  substituting  equations  (4-15),  (4-16)  and  (4-17)  into  the  model  equation 
(3-2)  , 

(2  e2  \ 

1  *  ~  -J  )  6in  6  +  K2  sin20  -  K3  sin3e 

+  (6  -  K  a  -  K  e  +  K.e  e.J  cos  0  -  £K  sin  0  cos  6  .  (4-18) 

'■l  1  12  1  3  4'  sp 


Comparing  equation  (4-18)  with  (3-12),  it  can  be  seen  that  the  measured  bias 

differs  from  the  true  bias  by  an  amount  K.e.e,  ,  the  difference  in  the  scale 

I  21  2 1\  114’ 

factors  is  2  +  e^/2j  and  the  axis  errors  differ  by  Kj (e2  -  •  The 

first  two  quantities  should  be  very  small,  depending  as  they  do  on  products  of 

-.3 

small  angles.  If  we  assume  Ej  ~  e^  ~  10  arc  minutes  (3  x  10  radians)  then  the 

bias  will  be  in  error  by  about  9  ug  and  the  scale  factor  by  about  9  ppm.  These 

may  be  considered  as  acceptably  small  errors  in  normal  testing;  furthermore,  e. 

-4 

is  more  likely  to  be  less  than  2  arc  minutes  (6  x  ]0  radians),  reducing  the 

bias  and  scale  factor  errors  to  <  2  pg  and  <  5  ppm  respectively.  The  dominant 
contributor  to  the  axis  error  uncertainty  is  the  KjC2  term  which  shows  that 
great  care  must  be  taken  to  level  the  mounting  surface  in  the  plane  of  rotation. 

An  error  e^  =  2  arc  seconds  is  equivalent  to  an  axis  error  change  of  10  pg/g  . 

A  similar  analysis  holds  for  mounting  position  2  in  which  the  pendulum  axis 
is  kept  nominally  parallel  to  the  axis  of  rotation.  Again  the  dominant  error 
angle  is  e2  which  determines  the  accuracy  with  which  the  appropriate  axis  error 
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is  measured.  Provided  that  the  other  dividing  head  and  mounting  fixture  errors 
are  kept  below  a  few  arc  minutes,  they  do  not  affect  the  precision  with  which 
accelerometer  parameters  are  determined. 

The  same  argument,  however,  is  not  true  for  the  third  orthogonal  mounting 
position;  that  in  which  the  accelerometer  is  rotated  about  its  sensitive  ai.is. 
In  this  case,  the  acceleration  components  along  the  true  axes  are 


a  =  (a  -  e.)  sin  0  +  (6  -  e.)  cos  6  -  e, 
s  4  3  1 


(4-19) 


/  -2  2  v 

■-B-*  * 


sin  6  -  (y  +  e0  ~  t.e.)  "os  0  +  e,e, 
2  3  4  1  4 


(4-20) 


2  2 


a  =  -  (v  +  e„)  sin  0  + 

n  2 


B-») 


cos  0  -  e^e^ 


(4-21) 


and  the  model  equation  (3-2)  becomes 


1  =  Ko~Ki£- 


]  ~  [O5,  ~  K,«)  +  V4]  Sin  6  +  +  K1 ^  "  Kle3] 


cos  0  .  (4-22) 


The  difference  between  the  true  and  measured  bias  is  Kj£  which  is  likely  to  be 
significant.  For  example,  if  £j  ~  5  arc  minutes,  then  the  bias  error  is 
~  1500  yg,  which  in  a  good  accelerometer  is  much  larger  than  itself!  This, 

however,  suggests  a  method  for  levelling  the  axis  of  rotation,  for  if  the  true 
bias  is  known  from  tests  in  one  of  the  other  mounting  positions,  then  the  axis 
orientation  can  be  adjusted  until  the  two  values  agree,  and  hence  =  0  . 

Measurements  of  the  components  of  axis  error  are  also  unlikely  to  agree 

with  those  obtained  in  the  other  mounting  positions  because  of  the  terms  K  e 

1  4 

and  K j £^  ,  both  of  which  can  be  large.  It  is  not  unusual  for  to  be  of  the 

order  of  10  arc  minutes,  thus  contributing  an  additional  3000  yg/g  to  the 
component  of  axis  error.  This  is  a  relatively  enormous  quantity  and  it  is 
difficult  to  envisage  how  £^  and  e^  could  be  reduced  to  the  few  arc  seconds 
necessary  to  make  measurements  in  this  mounting  position  meaningful. 

5  CONCLUSIONS 

In  this  Report,  a  mathematical  modelling  approach  has  been  used  to  explore 
methods  for  calibrating  precision  accelerometers.  While  the  :onventional  four- 
point  test  remains  an  important  means  of  monitoring  instrument  performance. 
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multiposition  tests  are  considered  superior  when  measurements  of  the  highest 
accuracy  are  required.  Because  such  tests  allow  the  accelerometer  characteristics 
to  be  examined  in  much  finer  detail,  it  is  possible  to  isolate  any  anomalous 
effects  and  to  decide  whether  the  ouosen  model  equation  is  appropriate.  If 
necessary,  different  models  may  be  tested  to  determine  the  best  fit  to  the  experi¬ 
mental  data.  An  additional  advantage  of  multiposition  testing  is  that  bad  data 
points,  caused  by  operator  mistakes  or  spurious  signals,  are  immediately  apparent 
from  the  plot  of  residuals. 

With  the  present  emphajis  on  measurement  precision,  errors  associated  with 
the  test  facility  must  be  carefully  considered.  Assuming  the  accelerometer  to  be 
mounted  on  the  dividing  head  so  that  it  rotates  about  one  of  its  non-sensitive 
axes,  then  it  has  been  shown  that  only  levelling  in  the  plane  of  rotation  is  of 
prime  importance.  Other  angular  errors  of  the  order  of  a  few  arc  minutes  can  be 
tolerated  without  significint  loss  of  accuracy.  However,  when  the  accelerometer 
is  rotated  about  its  sensitive  axis,  then  these  errors  have  a  profound  effect  on 
the  results. 
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